SIMULATIONS OF CORE COLLAPSE SUPERNOVAE 
IN ONE AND TWO DIMENSIONS 
USING MULTIGROUP NEUTRINO TRANSPORT 



A. MEZZACAPPA 
ORNL, Physics Division, Building 6010, MS 6354, P-0. Box 2008 
Oak Ridge, TN 37831-6354 

In one dimension, we present results from comparisons of stationary state multi- 
group flux-limited diffusion and Boltzmann neutrino transport, focusing on quan- 
tities central to the postbounce shock reheating. In two dimensions, we present 
results from simulations that couple one-dimensional multigroup flux-limited dif- 
fusion to two-dimensional (PPM) hydrodynamics. 



1 Introduction 

Current supernova modeling revolves around the idea that the stalled super- 
nova shock is reenergized by absorption of electron neutrinos and antineu- 
trinos on the shock-liberated nucleons behind itlil Key to this process is the 
neutrino transport in the critical semitransparent region encompassing the 
neutrinospheres. Potentiallyp|aiding this process is convection below thcj|icu.- 
trinospheres and the shockp-d^QEi although the verdict has been mixedH'Q'ErO'El 

In semitransparent regions, neutrino transport approximations such as 
multigroup flux-limited diffusion (MGFLD) break down. Moreover, neutrino 
shock reheating/revival depends sensitively on the emergent neutrinosphere 
luminosities and spectra, and on the neutrino inverse flux factors between the 
neutrinospheres and the shock. At the very least, computation of these quan- 
tities by approximate transport methods have to be checked against exact 
methods, i.e., against Boltzmann neutrino transport. 

We present results from one- and two-dimensional supernova simulations. 
In one dimension, beginning from postbounce slices obtained from Bruenn's 
MGFLD simulations, which are subsequently thermally and hydrodynamically 
frozen, we compute stationary state neutrino distributions with MGFLD and 
Boltzmann neutrino transport. From these, we compute and compare the 
neutrino luminosities, RMS energies, and inverse flux factors, and the net neu- 
trino heating rate, between the neutrinospheres and the shock. In two dimen- 
sions, beginning with the same postbounce profiles, we couple one-dimensional 
MGFLD to two-dimensional PPM hydrodynamics to investigate "prompt" con- 
vection below the neutrinospheres and "neutrino-driven" convection below the 
shock. At the expense of dimensionality, we implement neutrino transport 
that is multigroup and that computes with sufficient realism transport through 
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opaque, semitransparent, and transparent regions. 

2 Prompt Convection 

Without neutrino transport, prompt convection develops and dissipates in ~15 
ms in both our 15 and 25 Mq models. When transport is included, there is 
no significant convective transport of entropy or leptons in either model. Neu- 
trino transport locally equilibrates a convecting fluid element with its surround- 
ings, reducing the convection growth rate and asymptotic convection velocities 
by factors of 4-250 between p = 10^^~^^ g/cm^, respectively^ Our transport 
neglects the finite time it takes for neutrinos to traverse a convecting fluid 
element; consequently, our equilibration rates are too rapid, and transport's 
effect on prompt convection is overestimated. Nonetheless, it is difficult to see 
how prompt convection will lead to significant convective transport when these 
effects are taken into account. 

3 Neutrino-Driven Convection 

In our 15 Mq model, large-scale semiturbulent convection is evident below the 
shock and is most vigorous at i = 225 ms after bounce (we began our run 
at 106 ms after bounce when a well-developed gain region was present.) At 
this time, the maximum angle-averaged entropy is 13.5, and the angle-averaged 
radial convection velocities exceed 10^ cm/sec, becoming supersonic just below 
the shock. Despite this, at t = 506 ms after bounce our shock has receded, the 
convection below it has become more turbulent, and there is no evidence of an 
explosion or of a developing explosionQ 

Different outcomes obtained by the various groups are mainly attributable 
to differences in the neutrino transport approximations, which determine the 
postbounce initial models and the neutrino luminosities, RMS energies, and 
inverse flux factors, which in turn define the postshock neutrino heating rates. 
Most notable are dramatic differences in the RMS energiesQ For example, for 
Vi^D — r? ^^'^ T^j, =T{t^^ = 2/3), which is implemented by Burrows et al. 
(1995)0 the electron neutrino RMS energy and matter temperature at the 
neutrinosphere are related by < E'^^ >i/2— 3.6T, whereas our MGFLD cal- 
culations give < >i/2— 3.0T. Because the neutrino heating rates depend 
on the square of the RMS energies, the Burrows et al. rates would effectively 
be 40-50% higher. In the Herant et al. (1994)3 calculations, at 100 ms after 
bounce < E^^ >^ 13 to 14 MeV and < E^^ >~ 20 MeV, whereas we obtain 
the significantly lower values, 10 and 13 MeV, respectively. 
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4 Boltzmann Neutrino Transport 

Although differences in the MGFLD and Boltzmann transport luminosities and 
RMS energies between the neutrinospheres and the shock were seen, the inverse 
flux factors differed most. The fractional difference relative to MGFLD rose 
to ^ 20-25% just below the shock. The net heating rate showed an even 
more pronounced differences^ with the Boltzmann rate being 2-3 times higher 
just above the gain radiusQ (Small differences in the heating rates lead to 
large differences in the net heating rate near the gain radii, where heating and 
cooling balance.) These results have two ramifications: (1) It may be possi- 
ble to obtain explosions in spherical symmetry in the absence of convection. 
(2) Improved transport in two and three dimensions will give rise to more 
neutrino heating, and may give rise to more vigorous neutrino-driven convec- 
tion. Simulations in both spherical and axi- symmetry with one-dimensional 
Boltzmann transport are planned. We are also developing realistic multidi- 
mensional, multigroup transport. 
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